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Some recent results from several fixed-target dimuon production experiments at Fermilab are pre- 
sented. In particular, we discuss the use of Drell-Yan data to determine the flavor structure of the 
nucleon sea, as well as to deduce the energy-loss of partons traversing nuclear medium. Future 
dilepton experiments at RHIC could shed more light on the flavor asymmetry and possible charge- 
symmetry-violation of the nucleon sea. Clear evidence for scaling violation in the Drell-Yan process 
could also be revealed at RHIC. 



1. Introduction 

The first dilepton production experiment was performed 
at the AGS almost 30 years ago pj. Over the years, 
hadron-induced dilepton production experiments have 
led to the discoveries of various vector bosons (J/*, T, 
and Z°). They also provided important and often unique 
informations on parton distributions in nucleons, nuclei, 
and mesons. 

A series of fixed-target dimuon production experiments 
(E772, E789, E866) have been carried out at Fermilab in 
the last 10 years. Some highlights from these experiments 
are presented here. In addition, the prospect for perform- 
ing several dilepton production experiments at RHIC is 
discussed. 

2. Scaling Violation in Drell-Yan Process 

In the "Naive" Drell-Yan (DY) model, the differential 
cross section, m d 2 a / dx Fdm, has an expression showing 
its scaling property as follows: 
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T l a e lilo,(xi)q a (x2) + q a {x l )q a {x2)\. (1) 

The right-hand side of Eq. (1) is only a function of x\ , X2 
and is independent of the beam energies. This scaling 
property no longer holds when QCD corrections to the 
DY are taken into account. 

While logarithmic scaling violation is well established 
in Deep-Inelastic Scattering (DIS) experiments, it is not 
well confirmed in DY experiments at all. As an exam- 
ple, Figure 1 compares the NA3 data at 400 GeV 
with the E605 H and E772 @ data at 800 GeV. The 
solid curve in Figure 1 corresponds to NLO calculation 
for 800 GeV p + d (^/i = 38.9 GeV) and it describes the 
NA3/E605/E772 data well. No evidence for scaling vio- 
lation is seen. As discussed in a recent review [^|, there 
are mainly two reasons for this. First, unlike the DIS, 
the DY cross section is a convolution of two structure 
functions. Scaling violation implies that the structure 
functions rise for x < 0.1 and drop for x > 0.1 as Q 2 
increases. For proton-induced DY, one often involves a 



beam quark with x\ > 0.1 and a target antiquark with 
X2 < 0.1. Hence the effects of scaling violation are par- 
tially cancelled. Second, unlike the DIS, the DY experi- 
ment can only probe relatively large Q 2 , namely, Q 2 > 16 
GeV 2 for a mass cut of 4 GeV. This makes it more diffi- 
cult to observe the logarithmic variation of the structure 
functions in DY experiments. 

Possible indications of scaling violation in DY pro- 
cess have been reported in two pion-induced experiments, 
E326 §] at Fermilab and NAIO § at CERN. E326 col- 
laboration compared their 225 GeV 7r~ + W DY cross 
sections against calculations with and without scaling vi- 
olation. They observed better agreement when scaling 
violation is included. This analysis is subject to the un- 
certainties associated with the pion structure functions, 
as well as the nuclear effects of the W target. The NAIO 
collaboration measured 7r~ + W DY cross sections at 
three beam energies, namely, 140, 194, and 286 GeV. By 
checking the ratios of the cross sections at three different 
energies, NAIO largely avoids the uncertainty of the pion 
structure functions. However, the relatively small span 
in ^/s, together with the complication of nuclear effects, 
make the NAIO result less than conclusive. 

RHIC provides an interesting opportunity for unam- 
biguously establishing scaling violation in the DY pro- 
cess. Figure 1 shows the predictions for p + d at y/s — 
500 GeV. The scaling-violation accounts for a factor of 
two drop in the DY cross sections when ^fs is increased 
from 38.9 GeV to 500 GeV. It appears quite feasible to 
establish scaling violation in DY with future dilepton pro- 
duction experiments at RHIC. 

3. Flavor- Asymmetry and Charge-Symmetry- 
Violation of the Nucleon Sea 

The DY process complements DIS tool to probe 

parton distributions in nucleons and nuclei. This is well 
illustrated in the recent progress in the study of flavor- 
asymmetry in the nucleon sea. 

Until recently, it had been assumed that the distri- 
butions of u and d quarks in the proton were identical. 
While the equality of u and d is not required by any 
known symmetry, it is a plausible assumption for sea 
quarks generated by gluon splitting. As the masses of 
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FIG. 1. Comparison of DY cross section data with NLO 
calculations using MRST [15] structure functions. Note that 
r = X1X2. The E772 [4], E605 [3], and NA3 [2] data points 
are shown as circles, squares, and triangles, respectively. The 
solid curve corresponds to fixed-target p+d collision at 800 
GeV, while the dotted curve is for p+d collision at yfs = 500 
GeV. 



FIG. 2. The ratio of d/u in the proton as a function of x 
extracted from the Fermilab E866 cross section ratio. The 
curves are from various parton distributions. The error bars 
indicate statistical errors only. An additional systematic un- 
certainty of ±0.032 is not shown. Also shown is the result 
from NA51, plotted as an open box. The open circles cor- 
respond to d/u values extracted with the assumption of the 
CSV effect reported in ref. [10]. 



the up and down quarks are small compared to the con- 
finement scale, nearly equal number of up and down sea 
quarks should result. 

The assumption of u(x) = d(x) can be tested by mea- 
surements of the Gottfried integral S , defined as 
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where F$ and F£ are the proton and neutron structure 
functions measured in DIS experiments. The second step 
in Eq. |^ follows from the assumption of charge symmetry. 
Under the assumption of a symmetric sea, u p = d p , the 
Gottfried Sum Rule (GSR) §, I G = 1/3, is obtained. 

The most accurate test of the GSR was reported by the 
New Muon Collaboration (NMC) ||, which measured F% 
and F£ over the region 0.004 < x < 0.8. They deter- 
mined the Gottfried integral to be 0.235 ± 0.026, signif- 
icantly below 1/3. This result implies that the integral 
of d — u is nonzero. However, the x-dependence of d — u 
remained unspecified. 

The proton-induced DY process provides an indepen- 
dent means to probe the flavor asymmetry of the nucleon 
sea jlo| . An important advantage of the Drell-Yan pro- 
cess is that the x-dependence of d/u can be determined. 



The NA51 Q and the E866 |j collaborations have 
compared the proton-induced DY dimuon yields from hy- 
drogen and deuterium targets, and they deduced the ra- 
tios of d/u as shown in Figure 2. For x < 0.15, d/u 
increases linearly with x and is in good agreement with 
the CTEQ4M Q and MRS(R2) @ parameterizations. 
However, a distinct feature of the data, not seen in either 
parameterization, is the rapid decrease towards unity of 
d/u beyond x^ = 0.2. The E866 data clearly affect the 
current parameterization of the nucleon sea. The most 



recent structure functions of Martin et al. 15 1 (MRST) in- 
cluded the E866 data in its global fit and its parametriza- 
tion of d/u is shown in Figure 2. 

Many papers have considered virtual mesons as the ori- 
gin for the observed d/u asymmetry (see recent review of 
Kumano |16|). Here the ir + (du) cloud, dominant in the 
process, p — ► ir + n, leads to an excess of d sea. Compari- 
son of the E866 data with various theoretical models has 
also been made |i~7f . 

It should be emphasized that the extraction of d/u 
values from the NA51 and E866 DY experiments re- 
quired the assumption of charge symmetry, namely, d n — 
Up, u n = d p , etc. Evidence for a surprisingly large charge- 
symmetry-violation (CSV) effect was recently reported 
by Boros et al. @jn| based on an analysis of F2 struc- 
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FIG. 3. Predictions of a(p + p — > W + x)/a(p + p — > VF~:r) 
as a function of x_f at ^/s — 500 GeV. The dashed curve cor- 
responds to the d/u symmetric MRS SO' structure functions, 
while the solid and dotted curves are for the d/u asymmetric 
structure function MRST and MRS(R2), respectively. 

ture functions determined from muon and neutrino DIS 
experiments. A large asymmetry, d n (x) ~ 1.25u p (:r) for 
0.008 < x < 0.1, is apparently needed to bring the muon 
and neutrino DIS data into agreement. How would this 
finding, if confirmed by further studies, affect the E866 
analysis of the flavor asymmetry? First, CSV alone could 
not account for the E866 data. In fact, an even larger 
amount of flavor asymmetry is required to compensate 
the possible CSV effect . This is illustrated in Figure 
2, where the open circles correspond to the d/u values one 
would have obtained if the CSV effect reported by Boros 
et al. is assumed. Second, there has been no indication 
of CSV for x > 0.1. Thus the large d/u asymmetry from 
E866 for x > 0.1 is not affected. 

To disentangle the d/u asymmetry from the possible 
CSV effect, one could consider W boson production, a 
generalized DY process, in p + p collision at RHIC. An 
interesting quantity to be measured is the ratio of the 
p + p — ► W + +x and p + p — > W~ + x cross sections p0| . 
It can be shown that this ratio is very sensitive to d/u. 
An important feature of the W production asymmetry 
in p + p collision is that it is completely free from the 
assumption of charge symmetry. Figure 3 shows the pre- 
dictions for p + p collision at yfs = 500 GeV. The dashed 
curve corresponds to the d/u symmetric MRS SO' J2l]] 
structure functions, while the solid and dotted curves are 
for the d/u asymmetric structure function MRST and 
MRS(R2), respectively. Figure 3 clearly shows that W 



FIG. 4. Ratios of heavy-nucleus to deuterium integrated 
yields per nucleon for 800 GeV proton production of dimuons 
from the Drell-Yan process and from decays of the J/ip, ip', 
T(15), and T(2S + 3S) states [5]. The short dash and long 
dash curves represent the approximate nuclear dependences 
for the bb and cc states, A 96 and A ' 92 , respectively. 

asymmetry measurements at RHIC could provide an in- 
dependent determination of d/u. 

4. Nuclear Medium Effects of Dilepton Produc- 
tion 

From a high-statistics measurement of dilepton produc- 
tion in 800 GeV proton-nucleus interaction, the target- 
mass dependence of DY, J/^ff, and T productions 
have been determined in E772 [^2| As shown in Fig- 
ure 4, different nuclear dependences are observed for dif- 
ferent dilepton processes. While the DY process shows 
almost no nuclear dependence, pronounced nuclear ef- 
fects are seen for the production of heavy quarkonium 
states. E772 found that J/ "J and have similar nuclear 
dependence. The nuclear dependences for T, T' and T" 
are less than that observed for the J /if? and iff'. Within 
statistics, the various T resonances also have very similar 
nuclear dependences. 

Although the integrated DY yields in E772 show lit- 
tle nuclear dependence, it is instructive to examine the 
DY nuclear dependences on various kinematic variables. 
Using the simple A a expression to fit the DY nuclear 
dependence, the values of a are shown in Figure 5 as a 
function of xt{x2), M, xp, and pt- Several features are 
observed: 

1. A suppression of the DY yields from heavy nuclear 
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FIG. 5. Nuclear dependence coefficient a for 800 GeV p+A 
Drell-Yan process versus various kinematic variables [22]. 
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FIG. 7. The change of mean pj for nuclear target, 
A<p?) = (p 2 t)(A) - (p 2 t )(D), for 800 GeV p+A Drell-Yan pro- 
cess and J/\f r and T(IS') productions. The solid curves are 
best fits to the A-dependence of A(p|). The J/* and T(1S) 
productions have identical curves for the A-dependence fits. 
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FIG. 6. Nuclear dependence coefficient a for the Drell-Yan 
process [22] versus xf for (a) X2 < 0.075, right scale, and (b) 
X2 > 0.075, left scale. The thin solid lines show a = 1. The 
dashed line is a linear least-squares fit to the lower points. 
Also shown is the mean value of x\ for (b). 



targets is seen at small X2- This is consistent with 
the shadowing effect observed in DIS. In fact, E772 
provides the only experimental evidence for shad- 
owing in hadronic reactions. The reach of small xi 
in E772 is limited by the mass cut (M > 4 GeV) 
and by the relatively small center-of-mass energy 
(recall that x\X% = M 2 / s). p-A collisions at RHIC 
clearly offer the exciting opportunity to extend the 
study of shadowing to much smaller x. 

2. oi{xf) shows an interesting trend, namely it de- 
creases as xp increases. It is tempting to at- 
tribute this behavior to initial-state energy-loss ef- 
fect. However, there is a strong correlation between 
Xf and X2 (xf = x\ —x%)^ and it is essential to sep- 
arate the x f energy-loss effect from the x-i shadow- 
ing effect. Figure 6 shows a versus xf for two bins 
of X2, one in the shadowing region (x2 < 0.075) 
and one outside of it (x2 > 0.075). There is no dis- 
cernible xf dependence for a once one stays out- 
side of the shadowing region. Therefore, the ap- 
parent suppression at large xf in Figure 5 reflects 
the shadowing effect at small X2 rather than the 
energy-loss effect. 

3. cx(pt) shows an enhancement at large pt. This is 
reminiscent of the Cronin Effect |25| where the 
broadening in p t distribution is attributed to mul- 
tiple parton-nucleon scatterings. It is instructive 
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to compare the pt broadening for DY process and 
quarkonium production. Figure 7 shows A(p|), the 
difference of mean p\ between p-A and p-D inter- 
actions, as a function of A for DY, J/^>, and T(15) 
productions at 800 GeV. The DY and T data are 
from E772 while the J/V results are from 

E789 p7|, E771 @|, and preliminary E866 anal- 
ysis ]29||. More details on this analysis will be pre- 
sented elsewhere Figure 7 shows that (pi) is 
well described by the simple expression a+bA 1 ^ 3 . It 
also shows that the pt broadening for J/*f? is very 
similar to T, but significantly larger (by a factor 
of 5) than the DY. A factor of 9/4 could be at- 
tributed to the color factor of the initial gluon in 
the quarkonium production versus the quark in the 
DY process. The remaining difference could come 
from the final-state multiple scattering effect which 
is absent in the DY process. 

Baier et al. |30| have recently derived a relation- 
ship between the partonic energy-loss due to gluon 
bremsstrahlung and the mean p\ broadening accumu- 
lated via multiple parton-nucleon scattering: 

-dE/dz = | a s A(p 2 t ). (3) 

This non-intuitive result states that the total energy loss 
is proportional to square of the path length traversed by 
the incident partons. From Figure 7 and Eq. we deduce 
that the mean total energy loss, AE, for the p+W DY 
process is w 0.6 GeV. Such an energy-loss is too small 
to cause any discernible effect in the xf (or x\) nuclear 
dependence. As shown in Figure 6, the dashed curve 
corresponds to AE = 2.0 ± 1.7 GeV (for p+W), and 
the E772 data are consistent with Eq. || A much more 
sensitive test for Eq. || could be done at RHIC, where 
the energy-loss effect is expected to be much enhanced 
in A- A collision 
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